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ABSTRACT 

Diffuse 7-ray emission produced by the interaction of cosmic-ray particles 
with matter and radiation in the Galaxy can be used to probe the distribution 
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of cosmic rays and their sources in different regions of the Galaxy. With its large 
field of view and long observation time, the Milagro Gamma Ray Observatory is 
an ideal instrument for surveying large regions of the Northern Hemisphere sky 
and for detecting diffuse 7-ray emission at very high energies. Here, the spatial 
distribution and the flux of the diffuse 7-ray emission in the TeV energy range 
with a median energy of 15 TeV for Galactic longitudes between 30° and 110° 
and between 136° and 216° and for Galactic latitudes between -10° and 10° are 
determined. The measured fluxes are consistent with predictions of the GAL- 
PROP model everywhere except for the Cygnus region (/ G [65°, 85°]). For the 
Cygnus region, the flux is twice the predicted value. This excess can be explained 
by the presence of active cosmic ray sources accelerating hadrons which interact 
with the local dense interstellar medium and produce gamma rays through pion 
decay. 

Subject headings: gamma rays: observations 



Introduction 



One hundred years after their discovery, the origin, acceleration, and propagation of 
Galactic cosmic rays is still unclear. Supernova remnants (SNRs) and pulsars are the pre- 
ferred candidate sources of Galactic cosmic rays. The detection of TeV 7 rays and X-rays 
from the same locations within SN Rs provides strong evidence that electrons are accelerated 
in SNRs (lAharonian et al.ll2006al ). However, no compelling evidence for the acceleration of 
hadrons in SNRs has yet been found. The Galactic diffuse 7-ray emission originates from 
the interactions of cosmic rays (hadrons and electrons) with the matter and radiation flelds 
in the Galaxy. Cosmic ray hadrons interact with matter producing neutral pions which in 
turn decay into 7 rays while cosmic-ray electrons produce TeV 7 rays by inverse Compton 
(IC) scattering off the interstellar radiation flelds. Therefore, 7 rays can provide information 
about the density and spectra of cosmic rays throughout the Galaxy. Gamma rays above 
10 TeV, produced by the highest energy particles accelerated in SNR shocks, can be used to 
probe the acceleration processes in Galactic sources. The spatial distribution of these TeV 
7 rays can be compared to model predictions of the diffuse 7-ray production from tt'' decays 
and IC scattering and thus the relative contribution of hadronic and leptonic mechanisms 
can be investigated. 



EGRET data on the diffuse emission from the Galactic plane (iHunter et al.l 119971 ) 
show a harder 7-ray spectrum from the inner Galaxy t han predicted on th e basis of the 



cosmic-ray spectrum and intensity measured at Earth (IBertsch et al.l Il993l : IStrong et al. 
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200430) • Many explanations for this " GeV excess" have been proposed, including a harder 



Stroiig et al 



2004a 



cles ( Ide Boer et al 



Galactic proton or el ectron spectrum (IGralewicz et al.l 119971 : lAharonian &: Atoyaru 12000 



b;_ Porter &: Protherog 1997 ) and the annihilation of dark matter parti 



20051 ). Recent studies of the EGRET data have also investigated experi- 



mental uncertainties associated with the determination of the Galactic diffuse emission and 



are discussed in iMoskalenko et al.l (l2007al ) . IStecker et al.l (l200g ) concluded that t he EG RET 



Baughman et al.l (120071 ) con- 



sensitivity above 1 GeV may have been overestimated, while 
eluded it may have been underestimated, leaving the situation unclear. 

The first measurement of diffuse emission above 3.5 TeV from a large region of the 
Galacti c plane (Galactic longitudes 40° < / < 100°) indicated the existen ce of a TeV 



excess (lAtkins et al.l l2005l : iProdanovic et al.l l2007l : ICasanova fc Dingus! l2007l ). More re- 



cent me asurements of the diffuse emission near 12 TeV from the Cygnus region of the 
Galaxy (lAbdo et al.l l2007al ) also show an excess when compared to predictions of GAL- 
PROP, a numerical model of cosmic-ray propagation in the Galaxy. Recently, HESS has 
detected very high-energy (VHE) diffuse emission from the Galactic Center Ridge, that is 
correlated with giant molecular clouds. The spectrum of the diffuse emission from the Galac- 
tic Center Ridge is significantly harder thar i the spectrum of the d iffuse emission predicted 
by assuming the local cosmic-ray spectrum (lAharonian et al.ll2006bl ). These previous results 
from Milagro and HESS support the hypothesis that the cosmic-ray flux is likely to vary 
throughout the Galaxy. 



The Milagro (lAtkins et al.ll2004l ) experiment is a water-Cherenkov detector at an alti- 
tude of 2630 m. It is composed of a central 60 m x 80 m pond with a sparse 200 m x 200 m 
array of 175 "outrigger" tanks surrounding it. The pond is instrumented with two layers 
of photomultiplier tubes. The top, "air-shower", layer consists of 450 PMTs under 1.4 m 
of water while the bottom, "muon", layer has 273 PMTs located 6 m below the surface. 
The air-shower layer allows the accurate measurement of shower particle arrival times used 
for direction reconstruction and triggering. The greater depth of the muon layer is used to 
detect penetrating muons and hadrons. The outrigger array, added in 2003, improved the 
angular resolution of the detector from ~0.75° to ~0.45° by providing a longer lever arm 
with which to reconstruct events. Milagro's large field of view (~2 sr) and high duty cycle 
(> 90%) allow it to monitor the entire overhead sky continuously, making it well-suited to 
measuring diffuse emission. 

Here, the Milagro measurement of the diffuse emission around 15 TeV from a re- 
gion of the Galactic plane of longitudes / G [30°, 110°] and / G [136°, 216°], and latitudes 
b G [—10°, 10°] is presented. The measured 7-ray flux and the latitudinal and longitudi- 
nal profiles of the emission are reported and compared to predictions of the GALPROP 
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model Jstrong et al.lboool . l2004"al lbl: IPorter et aliboosh . In GALPROP, first the propagation 
of cosmic rays in tlie Galaxy is modeled, and tlien tlie 7-ray emissivities are calculated using 
the propagated spectra of cosmic rays and the gas and radiation densities. The "conven- 
tional" model is tuned to reproduce the local direct cosmic-ray measurements. The "opti- 
mized" model has been designed to reproduce the EGRET data by relaxing the restriction 
from the local cosmic-ray measurements. In this version of the model, the proton spec- 
trum is constrained by the cosmic-ray antiproton measurements and the electron spectrum 
is constrained using the EGRET data themselves. 

Below, the course of the analysis is described, followed by the presentation of the results 
and a comparison with GALPROP predictions. In the last section, likely interpretations of 
the observations are discussed. 



2. Analysis 



The Milagro data, c ollected between Jul y 2000 and November 2007, were analyzed using 
the method described in lAbdo et al.l (l2007al ). Only events with a zenith angle less than 45° 
are included, which corresponds to declinations between —7° and 81°. The event excess is 
calculated using the background estimation method described in I Atkins et al.l (120031 ) with the 
modification t hat the events are w eighted by a factor dependent on the 7-hadron separation 
parameter A4 ( lAbdo et al.ll2007al ). Only declinations < 70° are considered. This choice is 
governed by the fact that for 6 > 70° the Galactic equator turns parallel to the ri ght ascension 



axis. This causes the ratio of on- to off-time in the background calculation (lAtkins et al. 



20031 ) to become too big for signal bin sizes of 2 deg longitude by 4 deg latitude which is the 
bin size that is used in the Galactic longitude flux profile. As a result, the present analysis 
is insensitive for 5 > 70° or / G [111°, 135°]. 

Within the re gion studied here, Milagro has previously detected four sources and four 
source candidates ( Abdo et al. 2007a| jbl). The contribution from these sources is taken into 
account by modeling each source as a two-dimensional Gaussian plus a constant. The source 
location {RA,S), the amplitude and radial width of the Gaussian, and the constant are de- 
termined using a minimization. The excess from each source is then calculated bin by bin 
using the resulting Gaussian function and subtracted from the total excess in the 0.1° x 0.1° 
bin event excess map of the Galactic plane. The resulting diffuse event ex cess is conver t ed to a 
flux with a Monte Carlo simulation of extensive air showers (CORSIKA, Ifleck et al. I JlQQsh ) 
and of the Milagro detector (GEANT4. lAgostinelli et al. I (120031 )). The diffuse flux is calcu- 
lated assuming a power law photon spectrum with a differential spectral index a = —2.75. 
This spectral index was chosen to match the cosmic-ray spectrum in the energy range of this 
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analysis (around 10 TeV). For a spectral index of -2.75 the median energy of detected events 
used in this analysis is 15 TeV. 

Studies of possible sources for systematic errors have been performed. The size of 
the fit region around the eight sources and source candidates was varied. The Gaussian 
fits to the event excesses were performed in boxes centered around the sources of 4° x 4°, 
6° X 6°, and 8° x 8°. The flux determination was also repeated for spectral indices of 
-2.4 and -2.9. The variations of the calculated fluxes were found to be less than 18%. 
Another study concerned the background estimation. The background is estimated using 
the met hod of direct integra tion, calculating a background map using events over a 2hr time 
interval (lAtkins et al.ll2003l ). To account for possible contamination of the background from 
signal events two prominent regions in the Milagro sky map are excluded when doing this: 
a 2 by 2 degree region around the Crab Nebula location and a ±2.5 degree region around 
the Galactic plane. To investigate possible systematic effects in the background estimation 
the size of the region of exclusion around the Galactic plane was increased to ±5 degrees. 
The flux variations were found to be less than 7%. These systematic errors are added in 
quadrature to a 30% systematic error which is deriv ed from the under estimation of the 
Milagro trigger rate by cosmic ray shower simulations (lAbdo et al.ll2007al ). 



3. Results 
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Table 1. Gamma-ray emission from the Galactic plane around 15 TeV. 



Region Statistical Diffuse Flux ( x lO-^^^TeV-^cm-^s-isr-i) 

for |b| < 2" Significance GALPROP 

(1, deg) cr Milagro*^ optimized conventional 



30 < Z < 65 


5.1 


23.1 ±4.5;^:^ 


20.0 


4.9 


65 < 1 < 85 


8.6 


21.8 ± 2.51:^:8 


10.2 


2.7 


85<l< 110 


1.3 


< 7.1 (95% c.l.) 


5.8 


1.3 


136 < 1 < 216 


0.8 


< 5.7 (95% c.l.) 


3.1 


0.9 



^The first error represents the statistical, the second the systematic uncertainty. 
See text for details. 



- 7- 



The upper plot in Figure [T] shows the Galactic longitude profile of the 7-ray emission 
as measured by Milagro before and after source contributions are subtracted. A 7-ray flux 
increase towards the Galactic center is visible, as well as the Cygnus region {I G [65°, 85°]) 
with a "bump" in the flux profile, even after source contributions are subtracted. The source- 
subtracted data points in the lower plot of Figure [1] are overlaid with the 7-ray emission pro 



files as p redicted by the op timized GALPROP model, version 53_6102129RG (IStrong et al. 



)_p 

2004a| ]bl: iPorter et al.ll2008l ). The most significant discrepancy between model predictions 
and data appears in the Cygnus region. Table [T] lists the emission in four different regions 
with the Cygnus region being the most significant, 8.6 a above the background. The pre- 
dictions of the optimized and conventional GALPROP model are also given in Table [TJ 
The measured diffuse flux from the Cygnus region is two (eight) times higher than the op- 
timized (conventional) GALPROP prediction. The measured flux from the inner Galaxy 
(/ e [30°, 65°]) is consistent with the optimized model and about five times higher than 
predicted by the conventional model. Calculating the ratio of the flux measured in the inner 
Galaxy to the flux measured in the Cygnus region cancels systematic detector effects like the 
underestimation of the trigger rate. The flux ratio in the data is calculated to be 1.1±0.2 
(stat.) compared to a ratio of 2.0 predicted by GALPROP. For the regions / G [85°, 110°] 
and / G [136°, 216°], flux upper limits are quoted since the statistical significances above the 
background are less than two standard deviations. 

The energy spectra as predicted by the optimized GALPROP model in the Cygnus 
region and in the inner Galaxy, / G [30°, 65°], are shown in Figure [2] together with the 
EGRET and Milagro measurements. As can be seen, at Milagro energies the dominating 
GALPROP contribution to the diffuse flux is due to inverse Compton scattering of cosmic- 
ray electrons off the cosmic microwave background (CMB). The electron injection spectrum 
of the optimized GALPROP model is chosen such that the diffuse 7-ray emission spectrum 
matches the EGRET measurement at GeV energies. The injection spectrum is a broken 
power law with a spectral index of -1.5 below 20 GeV and a spectral index of -2.42 above 
20 GeV extending to a maximum electro n energy of 1000 T eV with an electron flux of 
1.4xlO-9MeV-icm-2s-isr^i at 34.5 GeV Jstrong et al.ll2004al jbll 



Figure [3] shows the Galactic latitude profiles for b G [—10°, 10°] in three regions between 
Galactic longitude I = 30 and I = 110. Gaussian fits (not shown) to the data distributions for 
/ G [30°, 65°] (left plot) and / G [65°, 85°] (middle plot) yield values for the mean consistent 
with 6 = 0° and for a of 0.9° ±0.3° and 2.0° ±0.2°, respectively. Assuming a larger exclusion 
region around the Galactic plane in the background estimation (see previous section) gives 
the same narrow width. The emission profiles as predicted by the optimized GALPROP 
model are overlaid. The blue line shows the total flux prediction, the green the inverse 
Compton, and the red the pion contribution. For both the GALPROP prediction as well 
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as the data the flux numbers become smaller and the distributions become wider farther 
away from the Galactic center. In order to compare the predicted profiles with the measured 
profiles, the iii both regions is calculated. In the inner Galaxy region (/ G [30°, 65°] 
and b e [—10°, 10°]) a of 18.3 with 20 degrees of freedom is derived, corresponding to 
a probability of 57% that the for a correct model exceeds the observed one by chance. 
In the Cygnus region, this chance probability is calculated to be 10"''. The discrepancy 
between the model prediction of the latitude profile and the data in the Cygnus region is 
investigated further by fitting the measured profile between b = 10° and b = —10° to a linear 
combination of the predicted pion and IC profiles, Cjc ■ Picib) + C^r ■ p-K(b). The factors C/c 
and C,r are varied independently between 0.1 and 10 in steps of 0.1. The minimum value 
of is obtained for C,r=6.9, i.e. an increase of the pion contribution with respect to the 
GALPROP prediction by a factor 6.9, and (7/(7=0.1. The resulting chance probability is 3%. 
Performing the same fit to the inner Galaxy latitude profile yields an increase of the pion 
contribution of C,r=5.8 and a decrease of the IC contribution of C7c=0.1 with respect to the 
GALPROP prediction. The chance probability of this result is 93%, showing that there is 
no significant improvement between the two cases. 

Figure [3] also shows the Galactic latitude profile for b G [—10°, 10°] in the region above 
Cygnus [l G [85°, 110°], right plot). No significant enhancement near the Galactic plane is 
visible. 



4. Discussion 

Measurements of the diffuse TeV 7-ray flux from the Galactic Plane as well as its spatial 
distribution, the latitude and longitude profiles, have been presented. The diffuse 7-ray flux 
was compared to predictions based on both the conventional and the optimized GALPROP 
models. In the / G [30°, 65°] range the optimized GALPROP prediction is consistent with the 
Milagro measurement. According to the GALPROP model the diffuse emission near 15 TeV 
is dominated by the inverse Compton component, which in turn is dominated by ~100 TeV 
electrons scattering off the CMB (60 to 70% of the IC component, see Figure [2]). Therefore, 
Milagro measurements can be interpreted as showing the first indication of the cosmic- 
ray spectrum up and beyond 100 TeV using the predictions of the GALPROP optimized 
model. The propagated average electron spectrum calculated by the optimized GALPROP 
model predicts a flux of about four times the locally measured flux and extends above 100 
TeV with a differential spectral index of -3. An alternative explanation would be a harder 
nucleon injection spectrum in the inner Galaxy than nieasured locally, but th i s would have 



to be checked against local antiproton measurements (IMoskalenko et al.lll998l : IStrong et al. 
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20071 ). Studies of the lateral 7-ray emission profile in the inner Galaxy show that the profile 
is rather narrow and suggest a better agreement with the pion assumption. Better sensitivity 
is needed however to significantly differentiate between the IC and the pion hypotheses. 

The Cygnus region is the region with the most significant fiux excess in the Milagro 
spatial maps. The emission from the Cygnus region of the Galaxy at longitudes I G [65°, 85°] 
shows an excess by a factor of two when compared to the GALPROP optimized model 
predictions. Also for the Cygnus region, the model predictions are dominated by IC scatter- 
ing of electrons off the CMB (see Figure [2]). However, the measured and predicted profiles 
are inconsistent (see Figure [3]). Decreasing the IC contribution (which has a broad profile 
distribution) and increasing the pion contribution (which has a narrow profile distribution) 
with respect to the model prediction improves the agreement in the shape, refiecting that 
the latitude distribution measured by Milagro is rather narrow compared to the distribution 
predicted by the optimized GALPROP model. 

The Cygnus region hosts intense star formation activity, and is abundant with molecular 
clouds and candidate cosmic ray sources. The HEGRA source TeV J2032+413 is located in 
the area with the largest concentration of molecular and atomic hydrog en in the Cygnus re - 
gion and spatially coincident with the Milagro source MGRO J203H-41 (lAbdo et al.ll2007al ). 
For TeV J2032-I-413 an association has been proposed with Cygnus 0B2, the largest clus- 



ter of more than 2700 identified young stars at 5000 light year distance (lAharonian et al. 



2005al : lAnchordoqui et al.l 120071 ). Very close t o TeV J2032+413 0, the VLA has detected 
a weak non-thermal shell supernova remnant (iButt et al.l 120061 ) and more recent XMM- 
Newton observations have revealed the presence of an extended X-ray source co-located 
with Te V J2032-I-4130, wh ich could be produced by an unknown population of faint X-ray 
sources (iHorns et al.l 120071) . In order to explain the X-ray and TeV emissions from TeV 
J2032+413, 



Horns et al 



(120071 ) proposed both a hadronic and a leptonic scenario. In the 
hadronic scenario, the X-ray emission would be produced by the synchrotron emission of 
secondary electrons and the emission above 10 TeV should show a hard spectrum. If the 
X-ray and TeV emissions were instead of leptonic origin, the spectrum above 10 TeV should 
be rather soft due to the unavoidable Klein-Nishina suppression of the inverse Compton cross 
section. 

If cosmic rays are injected into the Galactic interstellar medium by any of these candidate 
counterparts to TeV J2032-I-4130 and MGRO J2031+41, then the newly accelerated cosmic 
rays interact with the local gas and produce 7 rays within 100 pc from the source with spectra 
which might be significantly different from the avera ge 7-ray spectrum, because the cosmic- 
ray spectrum has not yet been steepened by diffusion (lAharonian fc Atoyanlll996l ). Assuming 
a distance of about 1 kpc, the extension of the entire Cygnus region is about 300 pc, and thus 
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a single accelerator might influence strongly the entire region. Assuming the diffusion process 
to be energy dependent, the emission from a molecular cloud close to such an accelerator 
could be due to high energy protons which have been accelerated first, have left the remnant 
first, and have already reached the cloud, as they diffuse faster. Thus this emission would 
have a flatter spectrum and could produce V HE 7 rays in the Milagro energy range during 



the fir st 10,000 years after the SN explosion (IGabici fc Aharoniaru 120071 : iMoskalenko et al. 



2007bl ). Leptonic mechanisms for the production of VHE 7 rays are disfavored by the fast 
cooling of such highly energetic electrons. Assuming a cosmic-ray source of total energy 10^^ 
ergs having exploded less than 10,000 years ago in the ISM and a mole cular cloud of total 
mass greater than 10^ and less than 10^ solar masses (jPobashi et al.lll994l ). the hard spectrum 
cosmic ray nucleons scattering off the targets in the Cygnus region would produce a flux 
between 10"^^ and 10"^^ TeV'^ cm ^s ^ near 15 TeV (see Figure 1 of iGabici fc AharonianI 



(I2OO7I )). Calculated from Table [T], the emission measured by Milagro near 15 TeV not 
explained by GALPROP is roughly 2.9 x 10"^"^ TeV~^cm~^s~^, meaning that only a few 
strong young accelerators in the Cygnus region are needed to explain the excess emission 
measured by Milagro. 

The results herein and previously presented by the Milagro Collaboration provide further 
evidence for the existence of cosmic-ray accelerators in the Cygnus region of the Galaxy , 
favoring hadroni c production mechanisms (iBeacom fc Kistlerll2007l : lAnchordoqui et al.ll2007l : 
Butt et al.l I2OO6I ). If the emission originates from hadronic interactions of hard spectrum 
cosmic rays, a corresponding flux of neutrinos will arise from decay of charged pions, and 
neutrino detect ors such as Icecube could provide the conclusive pro be of proton acceleration 
in the Galaxy (IBeacom fc Kistlerl 120071 : IGabici fc Aharoniaru 120071 ) . 



Experiments like GLAST (with its improved sensitivity and angular resolution with 
respect to EGRET) will be able to separate the truly diffuse 7-ray emission from a potentially 
unresolved source component, and probe the spatial distribution of the diffuse emission at 
GeV energies (the "GeV excess"). Experiments like the proposed High Altitude Water 
Cherenkov (HAWC) detector will be able to constantly survey large regions of the sky, in 
particular the Galactic plane, at 7-ray energies up to ~100 TeV with 10 to 15 times the 
sensitivity of Milagro. Because of its more southerly location it will also be more sensitive 
to Galactic plane regions below I = 65°. This will put stricter constraints on models like 
GALPROP and thus provide crucial information about the propagation of cosmic rays above 
100 TeV. 
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Fig. 1. — Galactic longitude profile of the 7-ray emission around 15 TeV in the Galactic plane 
as measured by Milagro. Upper plot: Red data points with dashed error bars - no subtraction 
of source contributions, black data points - after subtraction of source contributions. Lower 
plot: Source-subtracted profile overlaid with prediction of the optimized GALPROP model 
- the red line is the pion contribution, the green line the IC contribution, and the blue 
line represents the total flux prediction between Galactic latitudes ± 2 degrees. There are 
no data points in the region of longitude / G [—144°, 29°] because it is below the Milagro 
horizon. The region / G [111°, 135°] is excluded because the analysis method is insensitive 
here (see text for details). 
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Fig. 2. — Gamma-ray spectra of the diffuse emission as predicted by the optimized GAL- 
PROP model for the Galactic plane - left plot: inner Galaxy (/ G [30°, 65°]), right plot: 
Cygnus region (/ G [65°, 85°]). The red bars represent EGRET data, the black bar the 
Milagro measurement, where the length of the bar represents the statistical uncertainty 
only. The dark blue line represents the total diffuse flux predicted by the optimized GAL- 
PROP model, the dark gray line the extragalactic background, and the light blue line the 
bremsstrahlung component. The two contributions at Milagro energies are shown as red 
line, the pion contribution, and green line, the total IC contribution. The green dashed line 
shows the dominant IC contribution from scattering of electrons off the cosmic microwave 
background, which amounts to about 60 to 70% of the IC component at Milagro energies. 
Other IC contributions which are less important, such as infrared and optical, are not shown 
separately. 
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Fig. 3. — Source-subtracted Galactic latitude profile of the 7-ray emission around 15 TeV 
in the inner Galaxy (left plot), in the Cygnus region (middle plot), and in the region above 
Cygnus (right plot) as measured by Milagro (points with errors) and predicted by the opti- 
mized GALPROP model. The blue curve is the total 7-ray flux, the red curve the pion and 
the green curve the IC contribution. 



